Spectral Properties of Coupled Bose-Einstein Condensates 
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We investigate the energy spectrum structure of a system of two (identical) interacting bosonic 
wells occupied by N bosons within the Schwinger realization of the angular momentum. This picture 
enables us to recognize the symmetry properties of the system Hamiltonian H and to use them for 
characterizing the energy eigenstates. Also, it allows for the derivation of the single-boson picture 
which is shown to be the background picture naturally involved by the secular equation for H . After 
deriving the corresponding eigenvalue equation, we recast it in a recursive iV-dependent form which 
suggests a way to generate the level doublets (characterizing the H spectrum) via suitable inner 
parameters. Finally, we show how the presence of doublets in the spectrum allows to recover, in the 
classical limit, the symmetry breaking effect that characterizes the system classically. 
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I. INTRODUCTION 

After the recent observation (!]] of Bose-Einstein con- 
densation in dilute atomic gases realized by confining 
a macroscopic population of atoms in a potential trap 
|^| , a large amount of work has been devoted to con- 
struct devices where two condensates 0] are trapped in 
a double-well potential. The interaction of such coupled 
Bose-Einstein condensates (BEC) gives rise to quantum 
phenomena such as coherent tunneling and interference 
effects that have been the subject of a huge number of 
studies, both theoretic Q and experimental || . 

The present paper is focussed on the dynamical aspects 
inherent in the interaction of two condensates trapped in 
two identical wells. Such a dynamics has been studied 
thoroughly in Refs. H within the minimal interaction 
scheme 
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where the classical fields ipj(r, t) obey two coupled Gross- 
Pitaevskii equations (GPE), and U, v, T, describe the 
interatomic scattering, the external potential and the 
tunneling amplitude, respectively. Fields tpAr) (often 
called the wave functions of the condensate H) are de- 
fined as the expectation value ipj(r, t) = (^ (r, t)) of the 
field operators *5fj within the many-body quantum the- 
ory of BEC's 0. The negligible space dependence of 
tpj through the condensates allows one to set Ai/jj ~ 
so that Eqs. (^) reduce to a hamiltonian system with 
two complex degrees of freedom where the nonlinear cu- 
bic terms provide the system with an ample variety of 
interwell processes. 

Before discussing the goals of this paper, it is useful to 
briefly review the relevant traits of the dynamics issued 
from the space- independent form of Eqs. Due to the 
assumption that the bosonic wells are identical, the asso- 



ciated model Hamiltonian TL{ipx, 1P2) = C^(|V-'i| 4 + IV , 2| 4 ) — 
vAf-T^xtp^+^t) (A/"= l^iP + l^l 2 ) exhibits a per- 
mutational symmetry (PS) realized by the exchange of 
the dynamical variables ipi, ipz. Physically, this is repre- 
sented by the population exchange, |i/'j| 2 being the boson 
number of the j-th well up to a volume factor. 

The structure of the model phase space V (this is two- 
dimensional because M is a constant of motion) reflects 
in a nontrivial way the presence of the PS. For energies 
E < E*, (E* is a critical value of the energy depending 
on the model parameters) the orbits are placed concen- 
trically around the energy minimum, and cover a region 
Cq C V that has essentially the structure of the harmonic 
oscillator phase space. This can be shown to entail the 
oscillations of the two condensates' populations around 
the common value M j 2. Also, for E < E* the PS maps 
each orbit in itself. The remaining part of V , filled by 
orbits with E > E*, is formed by two (spatially) disjoint 
components C+, C_ such that V = C+ U C_ U Co. In 
this case any given energy value E is associated with two 
distinct orbits 7 + e C+, 7- € C_ such that j± — > 7 T 
under the PS action. The presence of two energy max- 
ima (located in C+ and C_ symmetrically) causes such 
a structure. The remarkable feature is that, when E 
crosses E„ from below, the system undergoes a symme- 
try breaking (SB) phenomenon |^| (governed by a bifur- 
cation mechanism ||) since the system 'must' choose to 
evolve either along 7 + or along 7_. Dynamically, this 
entails the emergence of a permanent gap between the 
condensates' populations. Such an effect is also called 
a self-trapping effect ||[l0) in that, within a finite range 
of E, the system never leaves the region C+ (C_) where 
was initially placed to go in C_ (C + ). 

The investigation at the quantum level of the scenario 
just described has been performed in Ref. |l0| through 
the model Hamiltonian 



H = U{n\ + nl) — vN — T(a ia + + a 2 a+) . 



(2) 



1 



which represents the quantum counterpart of the Hamil- 
tonian H(ipi,tp2) for the Laplacian-free Eqs. ([!]). A 
simple way to obtain H relies on the fact that for low 
numbers of bosons per well (a realizable experimental 
situation) one can replace the condensate wave functions 
ipj's with the raising (lowering) operator a*, (af) obeying 
the canonical commutators [<2;,ad~] = Sij, i = 1,2. Since 
[N,H] = 0, the total number of bosons N := n\ + n 2 , 
m := afdi, is a constant of motion (to simplify, we shall 
denote its eigenvalue by N as well) . The rigorous deriva- 
tion of H is effected in Ref. |l(J by using a mode depen- 
dent form of the field operator within the many-body 
quantum theory jll| of BEC's. It is important to recall 
as well that model (||) has been studied also in Ref. ||] 
from the viewpoint of dynamical system theory. 

In this paper we investigate the spectral properties of 
model (||) by combining the use of the PS and of a further 
symmetry involving the change T — ► —T. The latter will 
be called odd symmetry (OS) for recalling its basic role in 
determining the structure of the energy spectrum when 
the total boson number N is odd. Such symmetries are 
used extensively to show that: 

(?) after recovering the known nondegeneracy of the 
Hamiltonian spectrum (see Ref. [^|,^0)), each energy 
eigenstate is either symmetric or antisymmetric under 
the PS action, 

(m) the doublets (pairs of close energy levels occurring in 
the energy spectrum when the model parameters range in 
a suitable interval) always pair a symmetric eigenstates 
with an antisymmetric one; we shall show how this fea- 
ture plays a basic role in the classic limit, 
(Hi) the separation mechanism causing the splitting of 
the energy levels (the splitting effect has been observed 
numerically in Refs. |^,|l^]) can be explained in a purely 
analytic way, 

(iv) the fact that total boson number is even/odd dramat- 
ically influences the eigenvalues' parity under the OS. 

As to point (i), we wish to emphasize that the main 
consequence of the nondegeneracy is to prevent the SB 
phenomenon occurring in V as well as the ensuing self- 
trapping of the system on a specific orbit of the two ones 
associated with a given energy E > E*. This reflects the 
intrinsic tunneling effect due to the quantum nature of 
the system. 

For attaining results (i)-(iv) we first reformulate the 
initial quantum problem of two coupled (identical) wells 
through the Schwinger realization | 1^| of the spin alge- 
bra in terms of two-boson operators. This allows one to 
reconstruct Hamiltonian (|^) within its dynamical alge- 
bra (this is introduced in Appendix A). The form of H 
thus obtained can be interpreted in terms of a bosonic 
model defined on a nonhomogeneous linear lattice with 
one effective boson (single-boson picture). Such a picture 
has been derived in Ref. |l3| and is reviewed in Sec. II 
together with the related formal background. Thanks to 
its group-theoretic character such a formulation is appli- 



cable to many-well systems with any boson number N. 

In Sec. Ill the diagonalization of H is faced in a sys- 
tematic way by making explicit the PS and OS action on 
the components of the energy eigenstates. The resulting 
characterization of the eigenstates leads to identify the 
recursive expression of the eigenvalue equation of H both 
for even N and for odd N. This, in turn, enables us to 
recognize explicitly a hidden parameter able to control in 
an analytic way the level distance of the doublets. Such 
a parameter seems to suggest an alternative procedure 
to evaluate perturbatively the energy levels. A similar 
mathematical problem was analysed in Ref. B as to the 
problem of the dynamical tunneling through a separatrix, 
where the splittings of the doublets were traced by using 
the standard perturbation method. 



II. SPIN PICTURE OF THE TWO- WELL MODEL 

A convenient way to study the spectrum of H consists 
in reformulating H by means of the Schwinger picture of 
spin operators. The latter is a two-boson realization of 
the spin operators |L2] 
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satisfying the commutators [J a , Jb] = ie a bcJc of the alge- 
bra su(2), where a, 6, c = 1,2,3, and e a b c is the totally 
antisymmetric tensor [^2|. The generators J a of su(2) 
commute with the Casimir operator C = Jf + J| + Jf = 
J^(Ja + 1) which, in the Schwinger picture, leads to the 
identification J4 = (n 2 + n\)/2. Consistently, one can 
check that [J4, J a ] = 0. Therefore, su(2) can be used to 
rewrite Hamiltonian (0) which takes the nonlinear form 



H = 2[UJj - vJi + UJl - TJ X ] 



(3) 



The spin picture embodies explicitly in H the dimension 
(2J + 1) of the Hilbert space H(N) of physical states 
of H, where the eigenvalue J(= N/2) of J4 is the in- 
dex of the spin representation |l2[| . The standard basis 
B N = {\J,m),\m\ < J — N/2} (J 3 \J,m) = m\J,m)) of 
H(N) is related to the number operator states through 
\J,m) — \n\,n2), where J = (n 2 + n\)/2, m = (n 2 — 
n\)/2. Thus H and J a can be seen as (2J+ 1) x (2J+ 1) 
matrices. The fact that 2 J4 = TV is a constant of motion 
due to [H, N] = is ensured by [J 4 , J a ] = 0. 

An important consequence made explicit by the spin 
picture is that the nonlinear term J| arising in Eq. (j^) 
prevents the standard procedure of diagonalization of H 
via a unitary trasformation of the group SU(2). Such a 
procedure works only for matrices that can be written 
as linear combinations of the algebra generators in that, 
by construction, they can be reduced to one of the gen- 
erators by some appropriate group transformations [ fl2"| . 
Since the diagonalization is greatly simplified by such a 
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reduction, the latter practically identifies with the diag- 
onalization itself. When a matrix O contains nonlinear 
terms of its generating algebra (as in the case of Hamil- 
tonian (|J) with respect to su(2)), resorting to a larger 
algebraic structure A enables one to express O as a lin- 
ear combination of generators of A. Such an enlarged 
algebra is called a dynamical algebra Jl4j for O. 

About Hamiltonian H the problem is solved by rep- 
resenting the algebra su N (2) [subscript N{= 2 J) recalls 
the algebra link with the total boson number] as a sub- 
algebra of A = su(M) for a suitable values of M. In 
Appendix A we show that A = su^A+l). The realiza- 
tion in A of the su N (2) generators occurring in H reads 

J 3 = X q (J+l-q)n q , J+ =E 9 [q(2J+l-q)]H+b q+1 , 

J- = (J+) + (J± = Ji ± iJa), where 6+ (b q ) are rais- 
ing (lowering) bosonic operators, n q = b+b q , and q e 
[1, N + 1]. Recalling that the constraint Nf, := "E q n q = 1 
must be accounted for (see Appendix A) we shall call the 
realization just obtained single-boson picture. This fur- 
nishes the simplest way to represent in a linear form the 
nonlinear term of H . In fact, Hamiltonian (^) becomes 

H = C + 2im q [m 2 {q) n q - rR(q, J) (b+b q+1 + h.c.)} , (4) 

with r :=T/U, R(q, J) := [( J+l/2) 2 - (m(q) - 1/2) 2 ] 1/2 , 
m(q) := J + 1 — q, and C = 2[UJ 2 — vJ], where the 
quadratic term has been re-expressed as a linear combi- 
nation of number operators n q . A valuable alternative 
form of the original model (|J) is offered by Eq. (||) which 
recasts the two- well dynamics in terms of the dynamics of 
a single boson on a linear nonhomogeneous lattice. In the 
single-boson picture physical states are expressed as |l5| ] 

\9)=E 9 %b+\0,...,0) (5) 

with the normalisation condition Eq|-0 9 | 2 = 1 (see Ap- 
pendix A). The system dynamics thus takes place on a 
hypersphere inside C N+1 , where it can be interpreted in 
a classical form. This matches the fact that states (||) can 
be shown to be su^A+l) coherent states |l4|. Based on 
Eq. (||), the Schrodinger problem idt\^) = H\fy) can be 
rewritten as a set of equations of motion 

= 2Um 2 <£ m - T [R% +1 tf m+1 + R% * TO _i] , (6) 

where we have introduced ^> m := ip q , and Rm } :— R(q, J) 
with to :— J + 1 — q, to join the present formalism with 
the spin picture basis Sjv where states are labeled by 
rn G [—J, J]. It is worth noting how Eqs. @ can be 
derived as well from the effective Hamiltonian 

(H)=C +^ m [2Um 2 \^ m \ 2 -TR ( ^ ($* m * m _i+cc.)] , (7) 

representing the energy expectation value (\f'|ii"|\I f ) = 
(H), provided the Poisson structure {^m,^!} = S m e/ih 
is assumed (see commment pi). The time evolution of 



the dynamical variables (the components of | ) is deter- 
mined once the initial condition ^(O)) at t = has been 
assigned. Upon denoting by X m the components of the 
H eigenstates \X), one can retrieve the secular equation 

EX m — 2Um 2 X m — T \R^ +1 X m+1 + Rm X m -i\ , (8) 

from Eq. (Jfy. A m 's can be shown to be real. 

The procedure relying on the dynamical algebra con- 
struction has led to interpret model (g) as a lattice model 
with one boson via Eq. (^) . This procedure has been im- 
plemented as well for illustrative purposes since it shows 
clearly how the nonlinearity occurring in the matrix form 
of H is transferred to the coefficients of the linear com- 
bination in A. We point out that such a simplifica- 
tion also works for a linear chain of S interacting wells 
with TV bosons whose Hamiltonian can be written via 
the generators of su N (S). The latter, in fact, can be 
always immersed within an algebra su(M) with M suffi- 
ciently larger than 5*; moreover the two-boson realization 
of su(M) can be obtained for any M. 

The component form of secular equation (||) represents, 
at the operational level, a basic intermediate result. The 
latter, as shown in Sec. Ill, is used to characterize explic- 
itly (in the same spirit of Bloch's theorem for electronic 
wave functions) the inner spectrum structure as well as 
the structure of the energy eigenstates. 

III. SPECTRUM STRUCTURE 

In order to investigate the spectrum structure, we con- 
sider first the effects of the PS and the OS on the energy 
eigenstates, and make explicit how such two symmetries 
strongly characterize the eigenstate components. Then, 
we employ the results of such an analysis to recast the 
eigenvalue equation related to equation (Q) in a recur- 
sive form, and identify the parameters able to control the 
splitting of the energy levels. The case with half- integer 
J and integer J are treated separately. 

The PS is realized via the action of the unitary trans- 
formation U\ := exp[i7rJi] which takes J3 = (ri2 — ni)/2 
into UiJsU^ = — J3. This matches the effect of the PS 
classical action which involves the exchange of popula- 
tions n\ and n%. Let us introduce the hermitian operator 
P := exp[— in J] U\ whose action on the states \m) (we 
drop the representation index J in | J, m) since it is fixed) 
is deducible from the equation Ui\m) — exp[i7rJ]| — m). 
Owing to [H, Ui] = 0, P can be diagonalized together 
with H . The action of P on the standard basis, P\m) — 
I — to), implies that, for a generic state |\E'), 

P : I*) = S m * m |m) -» = E m * m | - to) . 

In particular, the P action on an eigenstate \X) entails 
the component transformation X m — > crA_ m , where a is 
not fixed, in Eqs. (pi). Actually, these remain unchanged 
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since R% +1 = R^ m for each to. The fact that P 2 = I 
fixes a showing how the allowed eigenvalues for P are 
a = ±1. Since P : X m — > ±X- m , each eigenstate has 
thus a definite symmetry character under to — > — to. 

This fact suggests to reorganize the vectors basis in 
terms of vectors |m, ±) = (|to) ± | — to))/\/2 that are 
either symmetric or antisymmetric. The new basis al- 
lows one to define in TL{N) a symmetric (antisymmetric) 
subspace H + (N) (H~(N)) whose vectors have compo- 
nents we will denote by <I> + ($~) such that P^ = ±$ ± . 
When the description in the new basis is adopted then 
Eqs. (||) for the eigenvector components can be written 
in the matrix form 



E 



x+ ' 




\ Sj(T) 





x- 







Aj(T) _ 



x+ 
x- 



(9) 



where Sj(T) (Aj(T)) is the sub-matrix associated with 
symmetric (antisymmetric) sector, and represents the 
zero-matrix. The matrix equation (JoJ) is separable in two 
independent equations for X + and X~: Their explicit 
form which depends on the representation index J is dis- 
played in the sequel. 

The odd symmetry (denoted by OS) is obtained by 
combining the action U3J1 = —J\ of U3 :— exp[z7rJ 3 ] 
on —TJ\ in H with the change T — > — T which restores 
the initial form of H. Since J can be either integer or 
half-integer, for considering the two cases separately it is 
convinicnt to introduce the matrix Lj(T) := ||L m f|| 

Lmi = 2Um 2 5 m ,i - T[R^S m ,i +1 + R ( e n 6 m +iA, (10) 

where m,£ = 1/2, 3/2, . . . , J if J is half-integer, and 
to, i = 1, 2, . . . , J if J is integer. 

Let us start with the half-integer case. In Eq. (||), the 
sub-matrices Sj (T) and A 1 (T) coincide with the matrix 
L j (T) up to the quantity —r/TR^L which must be added 
to the matrix element L11 with ?y = — 1(+1) in the an- 
tisymmetric (symmetric) case. Representing U3 in the 
basis ({|to, ±)}) entails 



U 3 






D~ 


-D 






(11) 



in which D — Diag(i, —i, . . .). The action of U3 on any 
vector takes its symmetric components into the antisym- 
metric ones and vice versa, namely PU3<& ± = ^U^^ 
if P^ = The structure of energy spectrum is 

reconstructed through the following three observations: 
i) The secular equation for X + and X~ derived from 
Eqs. (Q) can be written in the reduced form 

EC m = 2U m 2 C m — T \R { r ^ +1 C m +i + Rm C m -i] , (12) 

with C = X + ,X~, for 1/2 < to < J, whereas 

0= [2£/(l/2) 2 -vTR^L - E]C 1/2 -TR$ 2 C 3/2 , (13) 



holds for to = 1/2, where ij = — 1 (+1) in the antisym- 
metric (symmetric) case. For a given T / 0, Eqs. jl^ ) 
and ([[3]) show that the eigenvalues E can be seen as a set 
of functions E a (T,rj) of 77 with 1/2 < a < N/2, defined 
implicitly. The energy eigenvalue E a (T, +1) of each sym- 
metric eigenstate can be derived from that E a (T, —1) of 
an antisymmetric eigenstate by moving 77 from — 1 to + 1 , 
and vice versa. 

ii) We consider here the problem of ordering the set 
of eigenvalues E a (T, ±1). For T — > (this eliminates 
the 77 dependence) Eqs. (||) are solvable; the result- 
ing eigenvalues are doubly degenerate since the equa- 
tions for both the symmetric and the antisymmetric 
components are identical. Explicitly, T — » entails 
E a (T, +1), E a (T, — 1) — > 2Um 2 for some to which shows 
how the label a can be identified with m € [1/2, J]. 
The order induced by (positive) to on the set {2Um 2 : 
I mj < J} is inherited both by the symmetric eigen- 
values {E a (T, +1)} and by the antisymmetric eigenval- 
ues {E a (T,-l)} as proven by the limit T -> 0. This 
also implies that, for sufficiently small T, E a (T, ±1) 7^ 
E b (T,±l) iia^b. 

iii) Implementing the action of U3 whose matrix form is 
given by Eq. (O) on Eq. (Q) leads to the equation 



E 



x~~ 




\ Aj(-t) 





x+_ 







Sj(-T) 



x- 
x^ 



(14) 



where X^ = DX ± , and P X ± = ±X ± . The substitu- 
tion T — ► —T entails that [we use the simplified no- 
tation E±(T) := E a (T,±l)] the set of the symmetric 
(antisymmetric) eigenvalues {E+(T)} {{E~(T)}) coin- 
cides with the set of the antisymmetric (symmetric) ones 
{£ 6 -(-T)} ({££"(— T)}), where 1/2 < a, b < J. No- 
tice that, in general, E^(T) = E?(—T), where not nec- 
essarily b coincides with a. Nevertheless, for T — > 
E±{0) = E^(0) = 2Ua 2 = 2Ub 2 implies that b = a, 
as pointed out at point (ii). 

As a consequence of points (i)-(iii), we find that the 
symmetric eigenvalues are associated with the antisym- 
metric ones through the formula 



Et(T) = E*(-T) (1/2 < a < J) . 



(15) 



Also, since the eigenvalues equation can be cast in an 
iterative form via the recurrence equation 



d m (E) =(2Um 2 - E)d m+1 (E)-T 2 [Ri: n ^] 2 d 



+lJ Um+2 



(E), (16) 



which starts from 



= 



U/2-E + r,TR[ J / 2 d,(E) - T 2 [R[ J / 2 ) 2 d^E), (17) 



and terminates with dj(E) = 2UJ 2 — E, consistently 
with (iii) one finds that the eigenvalues cannot be even 
functions of T. 
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For integer J, the dimension of matrix Sj (T) changes 
from that of matrix Aj(T). In the antisymmetric case 
one finds Aj(T) = Lj(T), while in the symmetric 
case (5,7 (T) := ||S' m n||), where the indices runs over 
0, 1, . . . , J, one finds S i = S w = -TR { (\ S mn = L mn 
for m, n > 1, and S mn = otherwise. Eqs. ( |l2| ) still hold 
for integer J provided 2 < m < J. The two special cases 
are those corresponding to m = 0,1 







-ECq-TuR^ d, 



(18) 



= (2U - E) d - T [R ( 2 J) C 2 + <r R l ' } C ] , (19) 

with C = X + , X~. The parameter a must be set equal 
to one in the symmetric case (C = X + ), while in the 
antisymmetric case (C = X~) the expected elimina- 
tion of the component Xq follows from setting a = 0. 
Hence the dimensions of Hilbert sub-spaces are such that 
dim (N)= dim H + (N) — 1, while the secular equation 
for X~ will have a degree diminished of one. Explicitly, 
one has 



= 



E(2U-E) 2rn(J)l2 



d 2 (E)-E[R^] 2 d 3 (E). (20) 



In the symmetric case a ( J + l)-th degree equation for 
E issues from (||(]) through formula (0). Comparing the 
eigenvalue equations for the symmetric (er = 1) and anti- 
symmetric (er = 0) states shows that each, but one, sym- 
metric eigenvalue merges to an antisymmetric one when 
a goes from to +1. Due to the diversity of the secular 
equation with a = 1 from that with a = 0, even in the 
case with integer J, the energy spectrum is constituted 
by 2 J nondegenerate eigenvalues {E^(T) : 1 < a < J} 
that for T — > form J pairs Ef(T) — > 2Ua 2 and a single 
one Eq(T) which goes to zero in the same limit. Also, 
due to the quadratic dependence of Eqs. ([l6|), ( |20| ) on T, 
Eq. ( |l5| ) must be replaced with 



Et(T) = Et(-T) (l<a<J), 



(21) 



which, contrary to what happens with half-integer J, 
maps each eigenvalue in itself under T — > — T. In ad- 
dition, of course, one must consider Eq(T) = Eq(—T) as 
well. Figs. 1 illustrate the spectrum structure depen- 
dence on T/U for N = 6, 7 (see also Fig. 2). 



IV. DISCUSSION 

The interesting feature disclosed by the above analy- 
sis is the possibility to recognize both in the half-integer 
and in the integer case two inner parameters {rj and a) 
that control, in a way independent of T, the level split- 
ting generating the doublets. The limit T — > causes a 
coalescence of doublet levels such that E^(T) — » 2Ua 2 
which suggests T as a possible perturbative parameter for 



evaluating the level splitting. On the other hand, Fig. 2 
clearly shows that each eigenvalue E+(T) remains close 
to its partner E~ (T) on a finite range I a (T) of T indexed 
by the eigenvalue label a. In the half-integer case, this 
implies that, inside I a (T), indeed ij represents a good 
perturbative parameter (recall that E^(T) = E a (T,rf) 
with X] = ±1) which allows one to evaluate E a (T, ±1) by 
perturbing E a (T,r]), e. g., around rj — 0. For integer 
J, where the level separation is controlled by a G [0, 1], 
one can observe an effect similar to that showed in Fig. 
2: the symmetric eigenvalue E a (T, 1) = E+(T) remains 
close to its antisymmetric partner E a (T, 0) = E~ (T) on 
a finite range. Because the function E a (T,a) joins ana- 
litically E~ (T) to £J+ (T) then a can be assumed as the 
perturbative parameter for the present case. 

The actual size of the range I a (T) can be evinced 
roughly from Fig. 2, where the level separation strongly 
diverges only when Ef{T) cross E = E* := NT (recall 
that E* is the energy critical value defined in the intro- 
duction; its derivation can be found in Ref. fl3| ). This 
rule seems to be motivated from the insensitivity of Eq. 
([n]) (Eq. d20|)) from the parameter rj (a) in the terms 

r}TR[ J j 2 -E {a 2 [Rf] 2 -E 2 /T 2 ), 

for suitable values of E, T and of the other coefficients. 

Concerning the classical limit effected through N — > 
oo, numerical simulations with large boson numbers N 
show that the series of doublets becomes degenerate (co- 
alescence of the doublet levels) thus restoring the condi- 
tions that allow for the SB effect. How recovering the lat- 
ter is briefly illustrated via the following comparison be- 
tween the classical and the quantum behavior of the two- 
well system. Classically, at a given energy E > E*, the 
system described in V evolves either on 7_ or on 7+ (see 
Sec. I). Trajectories 7_ and 7+ are such that the pop- 
ulations' gap 7J2 — Tlx weakly oscillates around opposite 
values — fi and +fj,, respectively. Quantally, for E > E*, 
the combination |C^) = \X£) ± \X~) of the symmet- 
ric/antisymmetric eigenstates \X^) of the ath doublet 
can be shown to provide opposite nonvanishing expecta- 
tion values (J3) = ±x of J3 = (n 2 — n\)j2. This fact is 
caused by the eigenstate structure and was discussed in 
Ref. ||. Then, states \Cf) can be associated naturally 
to a pair of isoenergetic orbits 7_, 7+ that have fi = 2\- 
Increasing N, the time-dependent mixed state 

I*) = e- itE ZW n \X+) + e~ itE ^ T V n \X-) (22) 

(satisfying the Schrodingcr problem of H) exhibits a sort 
of temporary self-trapping effect (i. e. the localization 
either on 7_ or on 7+) which is repeated periodically and 
has a duration increasing with N . In fact, because of the 
oscillations of the factor exp{[it(E+ (T) - E~ (T)}} be- 
tween +1 and -1 entailing |*) oc |C+) and |\ff) cx \C~), 
respectively, the system stays in a quantum state involv- 
ing the localization on j± for intervals of the order of 
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the period At = 2iTh/[E+(T) - E~(T)] that increase 
when the level separation is reduced. The system re- 
mains definitively in the classical-like states (full emer- 
gence of the SB effect inducing the self-trapping) when 
the tunneling time from |C+) to \C~) diverges, namely 
for [E+(T) ~ E~(T)] (coalescence of doublet levels 
induced by N — > oo). 



V. CONCLUSIONS 

In Sec. I we have reviewed the dynamics of coupled 
Bose condensates described by Eqs. (|l|) (in the approx- 
imation with zero Laplacian terms) emphasizing the SB 
phenomenon that occurs in the phase space when increas- 
ing the energy over the critical value E*. Such a phe- 
nomenon (and the ensuing self-trapping effect on isoen- 
ergetic orbits j± £ V where ri2 — n\ oscillates around 
nonzero values ±//) has prompted the study of model (2) 
which represents the quantum counterpart of model ( 1 ) 
within the two-mode approximation of the condensate 
field operator. One of the purposes of the present work 
was to investigate the dynamical behavior corresponding, 
at the quantum level [through model (||)] , both to the SB 
effect and to the related self-trapping. An aspect we have 
particularly deepened is the quantum counterpart (level 
splitting) of the bifurcation mechanism generating pair 
of isoenergetic trajectories j± when E crosses E*. 

The formal set-up for studying the energy spectrum 
has been developed in Sec. II by recasting model (^) for 
the boson modes oi, a 2 into the matrix form (|3|) within 
the spin formulation a la Schwinger. The spin form of 
Hamiltonian (p|) makes easily viable the derivation of the 
secular equation Q . The latter has been achieved by us- 
ing the dynamical algebra method (this enacts systemat- 
ically the reduction of Hamiltonian nonlinearities) whose 
application is described in Appendix A. 

In addition to supplying equation (^), the use of the 
dynamical algebra method has shown the implicit link 
of the secular equation with the effective-bosons model 
(|J) . Such a model reformulates the two- well dynamics of 
Hamiltonian H with N bosons in a noticeably simplified 
form that consists of a single boson hopping on a non- 
homogenous lattice (single boson picture). The interest 
for the single boson picture and the underlying formal 
construction is motivated by the possibility of extending 
it to more structured models such as a chain model of S 
condensates with N boson. The case S — 3, which raises 
interest owing to its nonintegrable dynamical character, 
is presently under investigation |is| ] . 

Sec. Ill has been devoted to make explicit the struc- 
ture of the energy spectrum based on the symmetries of 
H . Upon introducing the pcrmutational symmetry (PS) 
and recognizing the further odd symmetry (OS), we have 
employed them to characterize the Hilbert space of H as 
well as the energy eigenstates, both for half-integer J 



(odd number of boson N) and for integer J (even num- 
ber of boson N). The A-dependent form of the secular 
equation obtained in the two cases has led to the central 
result of Sec. Ill, namely to recognize the possibility of 
introducing in a natural way inner parameters that con- 
trol the level splitting in the energy spectrum doublets. 

Such parameters -77 [a] is defined in Eq. ( p7| ) [Eq. (pf])] 
for half-integer [integer] J- have shown that the splitting 
originating the doublets can be generated explicitly in 
an analytic way. This fact combined with the doublets' 
structure exhibited in Fig. (2) suggests that parameters 
77, a can be used as perturbative parameters in approx- 
imating the doublet levels inside the regions of the T/U 
axis where the levels keep close. This approach may be 
preferable than the standard perturbative approach de- 
pending on the natural parameter T/U : this, in fact, can 
be shown to require higher and higher powers of T/U 
when approximating levels far from the ground-state |§ . 

We emphasize the fact that generating the level split- 
ting via the inner parameters 77, tr can be interpreted as 
the quantum form of the bifurcation mechanism issuing 
pairs of orbits j±. In general, our construction should 
furnish the quantum framework for describing the bi- 
furcation effects of any system whose Hamiltonian (in 
the critical regions of its phase space) has locally the 
same form of H. In this sense, both the quantum phase 
models and spin models JlTj in the mesoscopic system 
physics promise interesting applications. We notice as 
well that the matrix/ algebraic analysis underlying the 
study of 'quantum' bifurcation effects gives a valuable, 
both formal and practical, tool for characterizing quan- 
tally the chaos onset in the model with S = 3. 

The classical limit N — > 00 has been commented in Sec. 
IV, where the symmetry breaking (SB) effect inherent 
in the classical two-well dynamics is recovered from the 
quantum scenario via superposition ( ^) of the symmet- 
ric and antisymmetric states of each doublet. In fact, the 
coalescence of the doublet levels caused by N — > 00 (re- 
vealed by numerical simulations) leads, through a limit- 
ing process, to inhibit the oscillations of the system state 
I*) between |C+) and \C~) which ends up by coincid- 
ing with one of such states. This realizes the localization 
interpreted classically as the self-trapping effect. 

We conclude by illustrating a possible reformulation 
of Eq. (||) in a continuous form valid for large J di- 
rected to extimate the low part of energy spectrum. Set- 
ting X m =Y m /[(m + J)!(J -to)!] 1 / 2 in Eq. (f§ provides 
(2Um 2 -E)Y m = T[(J -m)Y m+1 + (J +m)Y m - 1 ], which 
reduces to the second order equation (F :— dF/dx) 



F - 2xsF + (2 J + 1 - sR + E/T)F/R = , (23) 

where s := ±1 and x := m/y/jR with R 2 := 1 + 2UJ/T, 
while Y(m) = exp[—am 2 /J] F(m/V JR), with a := 
(sR — l)/2, is the dependence of the (rescaled) compo- 
nents Y m on (the continuous variable) m. The assump- 
tion that F identifies with the n-th Hermite polynomial 
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gives the equation 2n = [2J+ 1 -R. + E/T]/R for s = +1 
which, in turn, supplies a set of eigenvalues. Their de- 
pendence on T is compared in Fig. 3 with the lowest part 
of the spectrum of N — 20 bosons. This result as well as 
the results/observations discussed above have prompted 
further work that is presently in progress. 



in a one-to-one correspondence with the states | J; to) , 
| m | < J = N/2, of the su N (2) standard basis (Ja\J; m) — 
m\J; to)) via the index map to = J + 1 — q. 
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APPENDIX A: 

The type of enlarged algebra we deal with is A = 
su(M) [the latter can be viewed as the generalized version 
of the spin algebra with M 2 — 1 generators] . Selecting 
an appropriate value of M allows one to rewrite the non- 
linear Hamiltonian H in terms of a linear combination 
of generators of A. This furnishes A with the status of 
dynamical algebra for H. To construct explicitly A it is 
useful to consider the two-boson form of the su(M) gen- 
erators E i;j := bfbj, (i ^ fhHij := [bjbi - bfbj)/2 that 
satisfy the commutators fl^j 

[Eij,Eki] — SjkEu — 5uEkj , 

[Eij,Hki] = {SjkEik — SjiEu + SnEij — 5ikEkj)/2 , 

with 1 < i,j,l,k < M. Within the present realization 
of su(M), the representation theory of semi-simple Lie 
groups |H states that the eigenvalue Q of the invariant 
operator Nb = Ej&^&j ([Nb,Eij] = 0) selects a specific 
representation of su(M) [Nf, can be viewed as the total 
particle number relatively to the creation (destruction) 
processes caused by bf (h)]. In fact, the dimension of 
the Hilbert space basis B(M,Q) = {\ni, tim)> Q = 
is given by 

dimB(M, Q) = {Q + M - 1)!/[(M - 1)\Q !] , 

where the states of the basis B(M, Q) are defined as 
\n±, ...,Um) = ®iL\\ n i) an d the number operator states 
\rij) fulfil the equations 

bi\rii) =^/rH\rii -1), bf\m) =Vni +T\th +1) • 

The simplest realization of A = su Q (M) is achieved by 
setting Q = 1 (single-boson picture) ; combining this fact 
with the requirement of preserving the dimension iV+1 of 
H(N), entails dimB(M, 1) = M = N + l thus providing 
A = sUj (N+l) as a dynamical algebra for H . The states 
of the related basis {\q) = 6+|0, ...,0),q £ [1,N + 1]} are 
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FIG. 1. First figure shows the energy eigenvalues vs. T/U for J = 7/2; second figure shows the case J = 3 (different types 
of line identify different eigenvalues) . 

FIG. 2. Solid lines represent the eigenvalues E^ vs. T/U for J = 15/2. The dashed line describes the curve E = 2JT. For 
each pair E^, the splitting becomes evident when 2JT reachs E^. 

FIG. 3. The figure shows the seven lowest ei gen values vs. T/U for J = 10 (solid lines); the seven dashed lines corrispond to 
the approximate eigenvalues derived from Eq. (G3) for s = +1. The approximation remains good for large values of T/U. 
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